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ABSTRACT 

We report a multisite photometric campaign for the (3 Cep stars V2052 Oph and 
V986 Oph. 670 hours of high-quahty differential photoelectric Stromgren, Johnson and 
Geneva time-series photometry were obtained with eight telescopes on five continents 
during 182 nights. Frequency analyses of the V2052 Oph data enabled the detection of 
three pulsation frequencies, the first harmonic of the strongest signal, and the rotation 
frequency with its first harmonic. Pulsational mode identification from analysing the 
colour amplitude ratios confirms the dominant mode as being radial, whereas the other 
two oscillations are most likely I — 4. Combining seismic constraints on the inclination 
of the rotation axis with published magnetic field analyses we conclude that the radial 
mode must be the fundamental. The rotational light modulation is in phase with 
published spectroscopic variability, and consistent with an oblique rotator for which 
both magnetic poles pass through the line of sight. The inclination of the rotation axis 
is 54° < i < 58° and the magnetic obliquity 58° < ^ < 66°. The possibility that V2052 
Oph has a magnetically confined wind is discussed. The photometric amplitudes of 
the single oscillation of V986 Oph are most consistent with an Z = 3 mode, but 
this identification is uncertain. Additional intrinsic, apparently temporally incoherent, 
light variations of V986 Oph are reported. Different interpretations thereof cannot 
be distinguished at this point, but this kind of variability appears to be present in 
many OB stars. The prospects of obtaining asteroseismic information for more rapidly 
rotating /3 Cep stars, which appear to prefer modes of higher Z, are briefly discussed. 

Key words: stars: variables: other - stars: early-type - stars: oscillations - stars: 
individual: V2052 Oph, V986 Oph - stars: magnetic field - stars: rotation 
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1 INTRODUCTION 

For over a century, the /S Cep stars are known to be vari- 
able on time scales of hours (Frost 1902), but it took half a 
century longer to understand the nature of their variability, 
radial and nonradial pulsations (Ledoux 1951). Nowadays, 
about 300 members of this class of pulsating star are known 
(Stankov & Handler 2005, Pigulski & Pojmanski 2008). 

Because of the simultaneous presence of radial and non- 
radial oscillation modes in these stars, and their rather sim- 
ple overall structure (basically a convective core and a radia- 
tive envelope), their potential as asteroseismic targets is evi- 
dent. Asteroseismology is the inference of the interior struc- 
ture of pulsating stars. This is accomplished by measuring 
their oscillation frequencies, comparing them with the eigen- 
frequencies of corresponding stellar models, and then fine- 
tuning those models to match the observed frequencies (see, 
e.g., Aerts et al. 2010, Handler 2012). 

Besides the Sun, the /9 Cep stars were the first main 
sequence pulsators for which clear constraints on their inner 
structure could be obtained asteroseismically (for a heavily 
abbreviated literature, see Aerts et al. 2003, Pamyatnykh 
et al. 2004, Handler et al. 2009, Aerts et al. 2011). Results 
indicate that further increases in heavy-element opacities are 
needed, and some stars have been shown to rotate faster in 
their interior than on the outside. 

These first successful studies were in most cases inten- 
tionally biased towards bright, slowly rotating stars. Slowly 
rotating /9 Cep stars driven by the k mechanism tend to have 
higher pulsation amplitudes (Stankov & Handler 2005) and 
therefore offer better possibilities for mode identification. 
Obviously, effects of rotation on the observed frequencies 
of axisymmetric modes of oscillation are also smaller, and 
rotationally split m-mode patterns would not overlap in fre- 
quency. This way of approaching asteroseismology of /? Cep 
stars proved to be sound. Therefore it appears reasonable to 
investigate targets that pose more difficult initial conditions, 
but that may also be more rewarding astrophysically. 

V2052 Oph (HR 6684, V = 5.8, B2IV-V) was discovered 
as a /3 Cep pulsator by Jerzykiewicz (1972), and its domi- 
nant mode identified as radial (Heynderickx, Waelkens & 
Smeyers 1994, Cugier, Dziembowski & Pamyatnykh 1994). 
Neiner et al. (2003) carried out an extensive multiwavelength 
spectroscopic and spoctropolarimetric study of V2052 Oph 
that revealed several interesting properties of this star. Be- 
sides the detection of a second, nonradial, pulsation mode, 
these authors could derive am accurate rotation period of 
3.638833 ± 0.000003 d. V2052 Oph also possesses a dipole 
magnetic field. Based on new data of superior quality, Neiner 
et al. (2012a) determined Bpoi « 400G, that the magnetic 
field is likely off-centred, and that He patches are present 
close to the magnetic poles. Because of the presence of a 
radial pulsation mode (that allows the determination of the 
mean stellar density), and of the known rotation period that 
makes it spin about twice as fast as the most "rapidly" ro- 
tating seismically well studied /3 Cep star (12 Lac, Desmet 
et al. 2009), it was deemed worthwhile to devote a large ob- 
servational effort to V2052 Oph. To this end, the present 
paper reports photometric results of a multisite campaign, 
whereas a companion paper (Briquet et al. 2012) deals with 
contemporaneous spectroscopy. 

Located only a few degrees in the sky from V2052 Oph is 



another (3 Cep star, V986 Oph (HR 6747, V = 6.1, BOHIn), 
with an interesting history in the literature. It is among 
the longest-period variables (P ~ 0.29 d, e.g., Jerzykiewicz 
1975) of its class, and among the most luminous and hence 
most massive (Jones & Shobbrook 1974, Stankov & Han- 
dler 2005). It is also a rapid rotator (wsini = 300kms~^ , 
Abt, Levato & Grosso 2002) and has been classified as a 
single-lined spectroscopic binary (Porb ~ 25.56 d, e = 0.23, 
FuUerton, Bolton & Penrod 1985). Frequency analyses pub- 
lished by different authors indicate variability with periods 
between 7 to 8 hours, but all studies noted that further 
photometric variability is present. However, no good expla- 
nation of its physical cause could be obtained (see Cuypers, 
Balona & Marang (1989) for a detailed discussion). Further- 
more, spectroscopic studies (FuUerton et al. 1985, Stateva, 
Niemczura & Iliev 2010) implied that the short period varia- 
tion is due to a mode of rather high spherical degree {1 — 4,6 
or 8). V986 Oph was also photometrically monitored during 
this multisite campaign, in the hope to gain understanding 
of its variability. 



2 OBSERVATIONS AND DATA REDUCTION 

Our photometric observations were carried out at seven dif- 
ferent observatories on five continents, from 11 March - 6 
September 2004. An overview of the campaign observations 
is given in Table 1. In most cases, single-channel differential 
photoelectric photometry was acquired but at Sierra Nevada 
Observatory a simultaneous uvby photometer was used. At 
observatories where no Stromgren vy filters were available 
we used Johnson V, with Stromgren u as possible comple- 
ment. Finally, as the photometer at the Mercator telescope 
has Geneva filters installed permanently, we used this filter 
system. In the three seasons preceding this campaign, 72.7 
hours of Geneva photometry had been obtained with the 
Mercator telescope. These are included here as well. 

We chose two comparison stars: HR 6689 {V = 5.96, 
A3V) was already used as a comparison for V2052 Oph by 
Jerzykiewicz (1972, 1993), but was in the second paper sus- 
pected to be variable. HR 6719 {V = 6.34, B2IV) was used 
as a comparison star by Jerzykiewicz (1993), but also sus- 
pected of variability. Subsequent photometric studies of the 
star did not mention variability, but Telting et al. (2007) 
reported line profile variations pointing towards pulsation 
of high azimuthal order. Although this choice of compari- 
son stars may not seem ideal, we did not find better suited 
candidates in this part of the sky. Fortunately, the com- 
parison stars proved to be constant within the accuracy of 
our measurements, and did not affect our results in any way. 
The targets were therefore observed alternatingly with these 
comparison stars, but V986 Oph was measured only in every 
other cycle. 

Data reduction was started by compensating for coin- 
cidence losses and subtracting sky background. Extinction 
corrections had to be made in two steps caused by the two 
comparison stars always being located at systematically dif- 
ferent air mass (similar right ascension, but different decli- 
nation and close to the celestial equator). This means that 
even small errors in the applied extinction coefficients cause 
variations in the nightly photometric zeropoints. 

Consequently, we first determined the extinction coefii- 
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Table 1 . Log of the photometric measurements of V986 Oph and V2052 Oph during the multisite campaign in 2004. Observatories are 

ordered according to geographical longitude. 



Observatory 


Telescope 


Amount 


of data 


Filter(s) 


Observer(s) 






Nights 


hours 






Tiibitak National Observatory, Turkey 


0.5m 


2 


6.1 


V 


TS 


South African Astronomical Observatory 


0.5m 


13 


50.3 


uvy 


EG, BN 


South African Astronomical Observatory 


0.75m 


7 


30.5 


uvy 


GH 


South African Astronomical Observatory 


0.5m 


9 


35.0 


uV 


TT 


Piszkesteto Observatory, Hungary 


0.5m 


3 


9.2 


V 


DL 


Sierra Nevada Observatory, Spain 


0.9m 


4 


19.7 


uvby 


ER 


Roquc dc los Mucliachos Observatory, Spain 


1.2m Mereator 


49 


155.2 


Geneva 


KU, MB, HW, GR, AM, CB, BV 


Fairborn Observatory, USA 


0.75m APT 


55 


198.7 


uvy 




Siding Spring Observatory, Australia 


0.6m 


40 


168.9 


uvy 


RRS 


Total 




182 


673.6 







cients with the standard Bouguer method from the compar- 
ison star measurements. We then examined the differential 
comparison star light curves for variability, resulting in a 
non-detection. Next, we imposed that the average nightly 
photometric zeropoints for each instrumental system be the 
same and correspondingly amended the extinction correc- 
tions within reasonable limits. This procedure considerably 
improved the accuracy of our final light curves, as examined 
with the target star data. The residual scatter in the differ- 
ential comparison star data is between 4.6 and 3.5 mniag in 
the Stromgren filters, and between 3.1 and 2.4 mmag in the 
Geneva measurements. 

Consequently, we computed differential light curves of 
the target stars and heliocentrically corrected their timings. 
The single-colour measurements were binned to sampling in- 
tervals similar to that of the multicolour measurements to 
avoid giving them higher weight in the consequent analyses. 
Finally, the photometric zeropoints of the different instru- 
ments were compared between the different sites and ad- 
justed if necessary. The resulting final combined time series, 
spanning 179.3 d, was subjected to frequency analysis. Light 
curves from the central part of our campaign are shown in 
Fig. 1, together with fits to be derived and commented on 
in what follows. 



Stromgren j// Johnson V light curves). The Stromgren b fil- 
ter measurements were not used because too few data are 
available. 

Amplitude spectra were computed, compared with the 
spectral window functions, and the frequencies of the intrin- 
sic and statistically significant peaJjs in the Fourier spectra 
were determined. Multifrequency fits with all detected sig- 
nals were calculated step by stop, the corresponding frequen- 
cies, amplitudes and phases were optimized and subtracted 
from the data before computing residual amplitude spectra, 
which were then examined in the same way. 

We consider an independent peak statistically signifi- 
cant if it exceeds an amplitude signal-to-noise ratio of 4 in 
the periodogram; combination signals must satisfy S/N > 
3.5 to be regarded as signiflcamt (see Breger et al. 1993, 
1999). The noise level was calculated as the average ampli- 
tude in a 5 d~^ interval centred on the frequency of interest. 

For the detection of pulsation frequencies we did not 
make use of the pre-campaign Geneva measurements as 
those were carried out without using comparison stars. They 
therefore have about a factor of 2.5 higher scatter than the 
campaign data and increase the noise level in a joint analy- 
sis. However, the pre-campaign measurements could in some 
cases be used to derive more precise frequency values. 



3 FREQUENCY ANALYSIS 

The heliocentrically corrected data were searched for peri- 
odicities using the program Period04 (Lenz & Breger 2005). 
This package applies single- frequency power spectrum anal- 
ysis and simultaneous multi-frequency sine-wave fitting. It 
also includes advanced options such as the calculation of 
optimal light-curve fits for multiperiodic signals including 
harmonic and combination frequencies. 

For purposes of frequency detection, the Stromgren u 
and Geneva U filter data were merged after checking that 
the oscillation amplitudes were the same within the errors. 
Measurements in the Stromgren y and .Johnson and Geneva 
V filters were treated as equivalent due to the same effec- 
tive wavelength of these filters, and were analysed together. 
After signals were believed to be detected, their presence 
was checked in the data of the individual filters (that is, 
the seven Geneva filters, Stromgren uv, and the combined 



3.1 V2052 Oph 

We started by computing the Fourier spectral window of 
the data, which turned out reasonably clean. The strongest 
aliases in the u/U and y/V data have only 36% of the am- 
plitude of the true signal. The amplitude spectrum of the 
data itself, dominated by the known radial mode frequency, 
is shown in the upper panel of Fig. 2. We chose to use the 
u/U data for presentation purposes for the pragmatical rea- 
son that all signals to be reported are detected in this data 
set alone. 

Prewhitening the strongest signal from the data and 
examining the residual amplitude spectrum, we recover the 
stellar rotation frequency (second panel of Fig. 2). Further 
analysis reveals two more signals in the frequency domain 
of the radial mode, as well as the first harmonic of the first 
mode and of the rotation frequency. The residual amplitude 
spectrum after prewhitening these six frequencies shows a 
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64.5 65 65.5 66 69.5 70 70.5 71 

HJD 2453100 + 



Figure 1. Some of our time-series photometry of V2052 Oph (upper three curves) and V986 Oph (lower three curves). The plus signs 
represent U or u measurements, the filled circles are v data and the open circles data in the V or y filter. The lines are the multifrequency 
fits derived from periodicity search. 



slight 1/f component, as expected from residual atmospheric 
effects in the data, and no signal in excess of 0.5 mmag. 

With frequency solutions for the individual filters as 
starting values, wc attempted to improve the accuracy of 
our frequency determinations by including the pre-campaign 
observations, therefore increasing the time base of the data 
set by a factor of 6.5. By examining the u/U and y/V data 
as well as confronting the results, we obtained more accurate 
values for all frequencies, with the exception of the second 
pulsation frequency where we encountered an aliasing prob- 
lem. Tests on which value would result in lower residuals etc. 
did not allow to determine a preferred value, and the choice 
of this frequency did not affect the outcome on the others. 
We therefore adopted the average of the two candidate val- 
ues, and used half the alias spacing as its uncertainty. The 
final values of the frequencies were then fitted to the cam- 
paign data alone but kept fixed, and only the amplitudes, 
phases and zeropoint were left as free parameters. The result 
of this procedure is listed in Table 2. 

The multifrequency fit listed in the table represents 
the data within rms residuals between 3.4 to 2.7 mmag 
(Stromgren data) and between 2.6 to 2.2 mmag (Geneva 
data). To search for possible additional signals, we merged 
the residual data from all filters and computed the combined 
amplitude spectrum (lowest panel of Fig. 2). It contains no 
peak in excess of 0.35 rrunag, and none with S/N ^ 3.2. 

Examining the wavelength dependence of the phases 
of the independent signals, we noticed that the dominant 
pulsation signal is not in phase in all filters, as demon- 
strated in Fig. 3. In particular, (f>v — (!>y = 0.9 ± 0.3° and 
<j)u — 4>y = 4.2 ± 0.2°, i.e. the shorter the wavelength, the 
later light maximum/minimum is reached. As the ampli- 
tudes of the other two signals in this frequency range are by 
at least a factor of 15 smaller, the errors in the phases axe 
correspondingly larger. Consequently, no statistically signif- 



icant phase shifts within the different filter passbands have 
been detected for /2 and /a. 

The Fourier parameters of the rotational light variation 
change substantially from filter to filter (cf. Table 2). To 
determine its shape wc first removed the pulsational vari- 
ability from the data, and then phased them with respect 
to the rotation period. The Geneva data were summed into 
20 phase bins and the more numerous Stromgren measure- 
ments in 25 bins. The rotational light curves are shown in 
Fig. 4. Unfortunately, these cannot be compared with coun- 
terparts from satelUte missions such as Kepler and CoRoT 
as our passbands are not sufficiently red sensitive. 

While we have arbitrarily phased these light curves and 
fits relative to H,JD 2453000.000, the hatched area in Fig. 
4 indicates the phase of minimum equivalent width of the 
ultraviolet spectral lines studied by Neiner et al. (2003). We 
will return to discuss this phasing in Sect. 5.1. 

3.2 V986 Oph 

The frequency analysis of our photometry of V986 Oph was 
performed in a similar way to that of V2052 Oph, and we 
also choose the u/U data for presentation. The amplitude 
spectrum of V986 Oph appears simple, with only one sig- 
nificant frequency present (Fig. 5). However, the residuals 
left behind a single-frequency solution arc between 6.3 and 
7.2 mmag per point in the Stromgren data, and between 7.0 
and 8.1 mmag per point in the Geneva data, thus about a 
factor of two to three higher than those for V2052 Oph or 
for the differential comparison star data. The poorer fit for 
V986 Oph is also readily visible in Fig. 1. The noise level in 
the residual amplitude spectrum of V986 Oph is even about 
a factor of five higher because of the smaller amount of data 
points available. Furthermore, the single frequency is not sig- 
nificantly detected in the Geneva data alone, neither in those 
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Table 2. Multifrcqucncy solution for our time- resolved photometry of V2052 Oph. Error estimates for the independent frequencies were 
derived from considering both formal errors (following Montgomery & O'Donoghue 1999) and differences in the u/U and y/V results. 
The quoted errors on the amplitudes are the formal values. The S/N ratio is for the u filter data. 



ID 


h 




h 




h 


2/1 


frot 


2 /rot 


Frequency (d^^ ) 


7.148474 ± 0.000005 


7.7567 


± 


0.0007 


6.82216 ± 0.00005 


14.296948 


0.27480 ± 0.00002 


0.54960 


u Ampl. (mmag) 


26.98 ± 0.10 


0.84 


± 


0.10 


0.54 ± 0.10 


0.52 ± 0.10 


2.76 ± 0.10 


0.78 ± 0.10 


V Ampl. (mmag) 


15.41 ± 0.09 


0.85 


± 


0.09 


0.60 ± 0.09 


0.20 ± 0.09 


1.97 ± 0.09 


1.05 ± 0.09 


y Ampl. (mmag) 


13.37 ± 0.08 


0.90 


± 


0.08 


0.56 ± 0.08 


0.23 ± 0.08 


1.62 ± 0.08 


0.29 ± 0.08 


U Ampl. (mmag) 


26.65 ± 0.16 


0.93 


± 


0.16 


0.97 ± 0.16 


0.53 ± 0.16 


2.89 ± 0.16 


0.47 ± 0.16 


Bi Ampl. (mmag) 


16.08 ± 0.15 


0.78 


± 


0.15 


0.97 ± 0.15 


0.26 ± 0.15 


2.56 ± 0.15 


1.35 ± 0.15 


B Ampl. (mmag) 


15.19 ± 0.15 


0.81 


± 


0.15 


0.75 ± 0.15 


0.27 ± 0.15 


2.21 ± 0.15 


0.91 ± 0.15 


B2 Ampl. (mmag) 


14.40 ± 0.16 


0.67 


± 


0.16 


0.60 ± 0.16 


0.07 ± 0.16 


1.88 ± 0.16 


0.34 ± 0.16 


Vi Ampl. (mmag) 


13.18 ± 0.15 


0.90 


± 


0.15 


0.76 ± 0.15 


0.38 ± 0.15 


1.47 ± 0.15 


0.29 ± 0.15 


V Ampl. (mmag) 


13.22 ± 0.13 


0.96 


± 


0.13 


0.74 ± 0.13 


0.29 ± 0.13 


1.40 ± 0.13 


0.32 ± 0.13 


G Ampl. (mmag) 


12.83 ± 0.16 


1.04 


± 


0.16 


0.60 ± 0.16 


0.41 ± 0.16 


1.56 ± 0.16 


0.27 ± 0.16 


S/N 


200.2 


6.4 


4.0 


3.9 


14.1 


4.0 



1 1 1 1 1 1 


: 

1 u data ~ 


^ rot u data - If 




2 u data - 2t -. 




Significance curve 



-| 1 1 1 1 — 

all data - 6f. 



i 



mmm 



5 10 

Frequency (d~^) 



15 



Figure 2. Amplitude spectra of our combined u/U filter data of 
V2052 Oph with consecutive prewhitening steps. The uppermost 
panel is practically identical with the spectral window of the data. 
The lowest panel shows the amplitude spectrum of the combined 
residuals in all filters. 
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Figure 3. Phase of the dominant oscillation mode of V2052 Oph 
with respect to that in the Stromgren y filter. The dots with 
the error bars are the measured phase shifts in all Geneva and 
Stromgren filters. 



obtained during the campaign, nor in the pre-campaign ob- 
servations. 

Investigating the matter deeper and keeping in mind 
that spectroscopic binarity of V986 Oph has been reported, 
we first looked for a possible light time effect. We there- 
fore merged the u/U, v and y/V data into bins no larger 
than 2.5 d as a compromise between not undersampling the 
reported 25.56 d orbit and having enough data points to 
determine the phase of the main light vaxiation. We did 
not find a statistically significant light time effect, within 
a limit of 1270 s at the orbital period. In one night of ob- 
servation (around HJD 2453158.9, see Fig. 1) a ~ 0.02 mag 
drop in light, suspicious of an eclipse, was present, but an- 
other data set obtained one prospective orbital period later 
did not show such a feature. 

Looking at the pulsation amplitude now, it appears that 
it dropped somewhat during the course of the campaign, but 
not exceeding the 2a level. We therefore assumed a constant 
amplitude for the remainder of this work, calculated the 
frequency of the single significant signal as a S/N- weighted 
average in the different Stromgren filters, and determined its 
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Figure 4. Phase diagrams of the rotational light variations of 
V2052 Oph in the different filters, plotted according to decreas- 
ing effective wavelength from top to bottom. Stromgren data are 
shown as open circles, Geneva data as filled circles, and fits are 
overlaid to guide the eye. The hatched area denotes the rotational 
phase zero as defined by Neiner et al. (2003). 
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Figure 5. Amplitude spectra and prewhitening of our combined 
u/U filter data of V986 Oph. 



Table 3. Frequency solution for our time-resolved Stromgren 
photometry of V986 Oph. The error on the frequency was de- 
termined from considering both formal errors (following Mont- 
gomery & O'Donoghue 1999) and differences between the u, v, 
and y frequency solutions. The quoted errors on the amplitudes 
are the formal values. The S/N ratio is for the u filter data. 





ID 


h 


Bl- 


Frequency (d~^ ) 


3.2886 ± 0.0003 




u Ampl. (mmag) 


4.0 ± 0.3 




V Ampl. (mmag) 


4.1 ± 0.3 




y Ampl. (mmag) 


3.3 ± 0.2 




S/N 


7.1 



amplitude and phase in all filters (see Table 3). The signal 
was found to be in phase within the errors in all passbands. 



4 MODE IDENTIFICATION 

We now attempt to identify the spherical degree I of the 
pulsation modes by means of the uvy and Geneva passband 
amplitudes of the pulsational signals detected in the light 
curves. Those amplitudes arc to bo compared with theoret- 
ically predicted ones from model computations, requiring 
the model parameter space to be constrained first. In other 
words, we need to determine the positions of the two target 
stars in the HR diagram as a starting point. 



4.1 V2052 Oph 

The latest spectroscopic Toff /log g values for V2052 Oph 
originate from Morel et al. (2006): Tcff = 23000±1000 K and 
log g = 4.0 ± 0.2, who also list v sin i = 61 kms^^ (including 
the macroturbulence velocity). Niemczura & Daszyiiska- 
Daszkiewicz (2005) derived T^s = 23300 ± 700 K and 
log g = 3.89 from low-resolution ultraviolet spectra. 

The online version of The General Catalogue of Pho- 
tometric Data (GCPD; Mermilliod, Mermilliod & Hauck 
1997) contains standard Stromgren and Geneva photomet- 
ric colours for the star. The Stromgren system calibrar 
tion by Napiwotzki, Schonberner & Wenske (1993), yields 
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Teff = 22700 ± 900 K, log 5 = 3.9 ± 0.3, and also provides an 
absolute magnitude estimate using the calibration of Balona 
& Shobbrook (1974): M„ = -2.59. The model atmosphere 
calibration of the Geneva system (Kiinzli et al. 1997) gives 
Teff = 22800 ± 500 K, log g = 3.8 ± 0.3. A relatively ac- 
curate HIPPARCOS parallax (van Leeuwen 2007) is also 
available: tt = 2.40 ± 0.41 mas. Adopting E{b - y) = 0.230 
from Stromgren photometry, this leads to Mv — — S.SIq ^. 

All these individual determinations are in very good 
agreement. We therefore assume Teff = 22900 ± 1000 K and 
log 5 = 4.0 ± 0.2. The tables by Flower (1997) then pro- 
vide BC = —2.20. A comparison of the Tofj/log g values 
with model evolutionary tracks prefers the lower value of 
our two absolute magnitude estimates, and suggests M = 
9.21°:? Me. 

Therefore we computed theoretical photometric ampli- 
tudes of the ^ I ^ 7 modes for models with masses between 
8.5 and 10.0 Mq in steps of 0.5 Mq, in a temperature range 
of 4.340 < log Teff < 4.379. We used OP opacities (e.g., 
Seaton 2005) and the Asplund et al. (2004) mixture. An 
overall metal abundance Z = 0.012 and a hydrogen abun- 
dance of X = 0.7 has been adopted, and no convective core 
overshooting was used. We are aware that Morel et al. (2006) 
and Niemczura & Daszynska-Daszkiewicz (2005) derived a 
somewhat lower metallicity, but this small inconsistency is 
not crucial in the mode identification process. 

We extracted theoretically calculated nonadiabatic pa- 
rameters from the models to determine the amplitudes in the 
different wavebands. This approach follows that by Balona 
& Evers (1999) and uses the same software, hence we re- 
fer to this paper for details on the procedure. Consequently, 
we computed the ratios of the theoretical amplitudes with 
respect to those in the Stromgren u filter, for modes of spher- 
ical degree ^5 / ^5 7 and frequencies between 6.3 and 8.3 
d~^ , and compared them with the observations (left-hand 
side panels of Fig. 6). 

The right-hand side of Fig. 6 shows a analysis of the 
photometric amplitudes, as outlined by Handler, Shobbrook 
& Mokgwetsi (2005) . These values use the measurements 
and standard errors of the amplitudes normalized to the 
mean of all passbands only, as the pulsation phases carry no 
additional information on the mode type in our case. 

Because of the high S/N of the dominant mode of 
V2052 Oph, the amplitude ratios in all the individual 
Stromgren and Geneva filters could be incorporated. Such an 
approach is not optimal for the two low-amplitude pulsation 
modes. For their identification, we considered the combined 
u/U, y/V as well as the Stromgren v data only. 

Like all previous authors, we identify the strongest 
mode as radial. The theoretically predicted amplitude ra- 
tios with respect to u are systematically higher than ob- 
served. We do not believe that this is a normalisation error 
because the measured u and U amplitudes agree quite well. 
We rather think that this is due to an imperfect choice of 
model parameters and because we ignored the chemical pe- 
culiarity of the atmosphere. Nevertheless, the data strongly 
suggest that the dominant mode is radial. 

The photometric amplitudes of the two other modes in- 
dicate they are of rather high and even spherical degree, 
most likely Z = 4 or 6. Our neglect of convective core over- 
shoot does not affect these mode identifications. Whereas 
overshooting would modify the deeper interior structure of 



the stax, the photometric amplitude ratios are mostly deter- 
mined in the stellar photosphere. 

Neiner et al. (2003) identified the 6.822 d"^ mode as ei- 
ther Z = 3 or 4. Therefore, the combined evidence points 
towards an Z = 4 mode, which is also the more likely iden- 
tification for the 7.757 d~^ signal because of reasons of ge- 
ometrical cancellation (Daszyriska-Daszkiewicz et al. 2002). 
This result is supported by Briquet et al. (2012), who could 
spectroscopically constrain the azimuthal order of the modes 
and who also discuss V2052 Oph more deeply in terms of 
convective core overshooting. 

4.2 V986 Oph 

For this star, no spectroscopic temperature and surface grav- 
ity determinations are available in the literature. The GCPD 
contains Stromgren colour indices for V986 Oph, but no 
Geneva indices. However, we have our own photometry in 
this system available. In Table 4, we first compare the stan- 
dard Geneva colours for V2052 Oph from the GCPD to 
those obtained from our data, demonstrating that they agree 
within a few millimagnitudes. Consequently, we trust the 
values that we obtained for V986 Oph. 

Using the standard values in the two photometric sys- 
tems as input for photometric calibrations, we must pro- 
ceed with caution. First, the Geneva colours of the star are 
somewhat out of the range of the calibrations by Kiinzli et 
al. (1997). Extrapolating their grids, one arrives at Teff w 
36000 K, log g ~ 3.8. Second, the mean Stromgren colours 
listed in the GCPD and the calibrations implemented by 
Napiwotzki et al. (1993) yield Tetr « 35600 K, log g « 3.0. 
The latter values however imply a stellar mass in excess 
of 4OM0, inconsistent with its BOIIIn spectral type. More 
reliable seems to be the Teff ~ 34700 K value from the cali- 
bration of the [u — b] index by Napiwotzki et al. (1993). 

Dciszyriska-Daszkiewicz (2001) determined Teff = 
30100 ± 2300 K, log 5 = 4.0 ± 0.5, [m/H] = 0.0, and 
E{B - V) = 0.228 ± 0.029 from lUE and visual fluxes. 
The value for reddening is consistent with E{b — j/) = 0.203 
from Stromgren photometry, but the effective temperature 
is much lower, and the surface gravity higher than from the 
photometric calibrations. These parameters rather imply a 
I6M0 star, but the large error bar on log g would allow 
masses up to 25Af0. 

For the purpose of example, we continue with Teff = 
34700 ± 1400 K, log 5 = 3.8 ± 0.3. We proceeded similar as 
we did for V2052 Oph, computing theoretical photometric 
amplitudes of the Z ^ 7 modes, but for models with 
masses between 25 and 36 Mq in steps of 1 Mq, and in a 
temperature range of 4.522 ^ log Tee ^ 4.558. A frequency 
range of 2.7 — 3.8 d~^ was considered for nonradial modes, 
and 3.2 — 3.4 d~^ for radial modes (to restrict the number 
of possible radial overtones). The comparison between the 
observed and computed amplitude ratios and a analysis 
are shown in Fig. 7, where we have restricted ourselves to 
the Stromgren data because the oscillation was not present 
at a significant level in the Geneva measurements. 

The results of this process clearly argue against a ra- 
dial pulsation mode. Considering nonradial modes, the ob- 
served amplitude ratios imply that the dominant signal in 
the light curve is most likely due to an Z = 3, 5 or 7 mode. 
The lowest x^i but also the smallest geometrical cancella- 
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Figure 6. Mode identifications for V2052 Oph from a comparison of observed and theoretical amplitudes in the Stromgren and Geneva 
bands. Left panels: amplitude ratios, normalised to unity at u. The filled circles with error bars are the observed amplitude ratios, the 
thin error bars denote the uncertainties in the theoretical amplitude ratios. The full lines are theoretical predictions for radial modes, 
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analysis of the photometric amplitudes. The smaller x^, the more likely an identification is. 



tion, then favour an identification as I — 3, if taken at face 
value. Adopting a lower mass, as implied by the Tcs/logg 
values by Daszyhska-Daszkiewicz (2001) results in a quali- 
tatively consistent picture with / = 3, 5 or 7 as the modes 
best reproducing the observed amplitude ratios. We refer to 
the discussion of the credibility of this mode identification 
near the end of Sect. 5.2. 

5 DISCUSSION 
5.1 V2052 Oph 

As mentioned in the Introduction, the presence of a radial 
mode in the star's pulsation spectrum allows to derive its 
mean density, provided the radial overtone is known. To 



this end, model evolutionary tracks were computed with 
the Warsaw-New Jersey stellar evolution code, for a rota- 
tional velocity of 80 kms~^ on the ZAMS (to match the rota- 
tion period) and other input parameters as specified in Sect. 
4.1. Nonadiabatic mode frequencies were calculated with the 
Warsaw pulsation code (e.g., sec Pamyatnykh et al. (1998) 
for a description of these codes), and models sought that 
had a radial mode at the observed frequency. Fig. 8 shows 
the result of this procedure in the form of a theoretical HR 
diagram. 

The effective temperature and luminosity derived for 
V2052 Oph in Sect. 4.1 is in best agreement with the hy- 
pothesis that the radial mode is the fundamental, but it 
cannot be excluded that it is the first overtone. In the first 
case, the stellar radius would be 5.3 ± 0.1 R© and the in- 
clination of the rotation axis thus 54° < i < 58°, given 
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Table 4. Geneva visual magnitudes and colours^ for our target stars. 



Star VM U V Bl B2 VI G 

V2052 Oph (literature) 5.803 0.598 0.855 0.834 1.543 1.563 2.018 

V2052 Oph (this work) 5.805 0.594 0.856 0.832 1.543 1.561 2.014 

V986 Oph (this work) 6.119 0.268 0.984 0.791 1.584 1.682 2.163 



^VM is the visual magnitude, whereas the other parameters are colour indices with 
respect to the B band magnitude (Golay 1972). 



the star's rotation period and v sin i. These values change to 
6.45 ± 0.15 Ro and 42° < i < 44° if the radial mode was 
the first overtone. Unfortunately, the scarcity of additional 
pulsation modes and their unknown azimuthal order leave 
meager prospects for asteroseismic constraints other than 
deriving the mean stellar density. 

Attempting to fit the I — 4 mode frequencies with mod- 
els of the same radial mode period but varying mass gave 
a number of possible solutions. Unsurprisingly, if the radial 
mode was the fundamental, models of lower mass appear 
more likely because these are more evolved (cf. Fig. 8) and 
therefore have more mixed modes of Z = 4. If the radial mode 
was assumed to be the first overtone, no such preference was 
seen. 

The only statement we can make is that the two Z = 4 
modes are unlikely to be rotationally split m-modes of the 
same radial overtone unless allowing for differential interior 
rotation. The same conclusion was reached by Briquet et al. 
(2012), with their independent model approach and identi- 
fications of m. The reason for this finding is that the effect 
of the Coriolis and centrifugal forces on the rotational fre- 
quency splitting arc very similar at the given rotation rate 
and for Z = 4, no matter whether the mode under consider- 
ation is a p, g or mixed mode. 

Neinor et al. (2003) determined a very precise rotation 
period for the star despite a non-optimal temporal distribu- 
tion of their data in terms of annual aliasing problems. The 
rotation period we have obtained is consistent within the er- 
rors with the more precise one by Neiner et al. (2003), and 
accurate enough to rule out that their value is affected by 
aliases. Therefore we confirm 3.638833 d as the best avail- 
able rotation period of V2052 Oph. 

In Fig. 4, the phase of minimum rotational light vari- 
ation coincides with minimum magnetic field strength and 
minimum UV spectral line equivalent width, as determined 
by Neiner ct al. (2012a). In Bl,v and B there is also a 
double maximum, as in the UV line strength. The shape of 
these variations indicates that both magnetic poles are seen 
during a rotation cycle, i.e. the sum of the angles of the 
inclination of the rotation axis and the magnetic obliquity 
i + l3 must exceed 90° . 

Neiner et al. (2012a) determined i to be 53° < i < 77° 
from modelling Stokes profiles. This is consistent with the 
inclination of the rotation axis we obtained with a stellar 
radius corresponding to radial fundamental mode pulsation, 
but not with the value assuming the radial mode is the first 
overtone. We therefore rule out the latter possibility. Using 
r = Bmin/Smax = cos(i — P) / cos{i + jS) , where r is the ratio 
of the minimum and maximum magnetic field strength of 



the rotation cycle, and the r values determined by Neiner et 
al. (2012a), we obtain 58° < ^ < 66°. 

As mentioned, the light curve of V2052 Oph shows rota- 
tional modulation (Fig. 4) in phase with the magnetic field 
and UV wind variations. This could be due to spots at the 
surface of the star (such as those suggested by Neiner et 
al. 2012a) or magnetically confined clouds in the circumstel- 
lar envirounic^iit (e.g. Townsend & Owocki 2005). Using the 
magnetic field value (Bpoi = 400 G, Neiner et al. 2012a), 
the wind velocity estimated in UV data [vint = 500kms~^ , 
Neiner et al. 2003), the stellar parameters from Sect. 4.1 
and above (M = 9.2 Mq, R = 5.3 Hq, vsini = 61kms^^ , 
i = 56°), and a mass loss typical of a B2 star (M = 10~® 
Mq yr~^), we derived the magnetic confinement parameter 
r?« (see ud-Doula & Owocki 2002) of V2052 Oph. 

We obtained that rj, = 2576, the Alfven radius is 
Ra = 7.12R, and the Kepler radius is Rk = 3.05i?. This 
implies that magnetic confinement should occur (r)* > 1) at 
the magnetic equator between Rk and Ra- Indeed in this 
region wind particles get trapped in closed field loops and re- 
mains centrifugally supported. Above Ra material escapes 
as the wind streches field lines open. Below Rk material 
lacks sufficient centrifugal support and falls back onto the 
star, but this transient material can still create a dynamical 
magnetosphere (ud-Doula, Owocki & Townsend 2008, Petit 
et al. 2011). 

However, a centrifugally supported magnetosphere usu- 
ally produces Ha emission and such emission has never been 
observed in V2052 Oph. Moreover, Oskinova et al. (2011) 
showed that V2052 Oph is only very weakly X-ray lumi- 
nous. A centrifugally supported magnetosphere could exist 
without producing Hq or much X-ray emission if the den- 
sity or temperature of the wind was not appropriate or if the 
confinement timescale was too long (see, e.g., Neiner et al. 
2012b for a more detailed discussion of the emission measure 
in magnetospheres). 

V2052 Oph is not the only magnetic (3 Cep star known. 
Telting, Aerts & Mathias (1997) spectroscopically detected 
rotationally equally split frequencies around the dominant 
radial pulsation mode of /3 Cep itself, and discussed whether 
these could be temperature spots on the surface or due 
to a magnetically distorted oblique pulsation mode (the 
magnetic field was reported by Henrichs et al. 2000). On 
the other hand, no such rotational frequency splitting has 
boon reported for CMa (Saoson, Briquet & Aerts 2006, 
Fourturie-Ravard ct al. 2011). Both stars have a single or 
dominant radial mode, such as V2052 Oph. For our target, 
we find no signals at frequencies split by one or two times 
the rotation frequency around the radial pulsation, within a 
limit of 0.2 mmag in amplitude. 
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Figure 7. Mode identifications for V986 Opli from a comparison 
of observed and tficorctical uvy amplitude ratios, normalised at 
u. In the upper two panels, the filled circles with error bars are 
the observed amplitude ratios. The full line in the uppermost 
panel is the theoretical predictions for radial modes, the dashed 
line for dipole modes, and the dashed-dotted lines for quadrupole 
modes. In the middle panel, the full line is for I = 3 modes, 
the dashed line for I = 4, the dashed-dotted lines for I = 5, the 
dotted lines for I = 6, and the dashed-triple-dotted lines are for 
I = 7 modes. The thin error bars denote the uncertainties in the 
theoretical amplitude ratios. The lowest panel shows the results 
of a analysis of the amplitudes. 



5.2 V986 Oph 

In Sect. 3.2 wc remarked that the residual scatter in our 
light curves prewhitened by the single coherent frequency is 
considerably higher than that in the other time series from 
this campaign. Since the two comparison stars are farther 
apart from each other on the sky as V986 Oph is from either 
of them, this cannot be due to residual data reduction er- 
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Figure 8. Constraints on the position of V2052 Oph (star sym- 
bol with error bars) in the HR Diagram. Some stellar evolutionary 
tracks are plotted and labelled with corresponding masses, and 
the theoretical edges of the /3 Cep (dashed-dotted line) and Slowly 
Pulsating B (SPB) star (dotted line) instability strips (Pamyat- 
nykh & Ziomek 2007) are shown. The thin full lines connect mod- 
els with radial modes at the observed frequency; "F" stands for 
the fundamental mode, " IH" for the first overtone. 



rors. Furthermore, the scatter in colour light curves of V986 
Oph, e.g. M — y, is much loss than the residual scatter in the 
individual passbands, and virtually the same as in the differ- 
ential comparison star data in the same filter combination. 
Wo therefore conclude that this high apparent scatter in the 
light curves actually represents intrinsic variability of V986 
Oph, and that this kind of variability is not dominated by 
changes in the stellar effective temperature. 

The main variability frequency we found is consistent 
with the one in the 1987 data by Cuypers et al. (1989), and 
the amplitude is comparable. However, these authors also 
remarked on the presence of longor-torm light variations on 
time scales longer than half a day. We do not find coherent 
variability on such time scales in the complete data set. We 
therefore subdivided our measurements into chunks compa- 
rable to the extent of the data by Cuypers et al. (1989) and 
analysed the residuals after prowliitening the main period- 
icity. Longer-term variability is present, but we find nothing 
periodic. Unfortunately, our data have too low a duty cy- 
cle for invoking techniques such as Time-Fourier analysis to 
search for some possible short-lived periodic variations. 

V986 Oph is not the first case of a /? Cep star for which 
unknown additional variability besides that of pulsational 
origin heis been found. Jerzykiewicz (1978) and Handler et 
al. (2006) discussed this problem for the star 12 (DD) Lac 
(11.5 Mq), Jerzykiewicz et al. (2005) for v Eridani (9.6 M©), 
and Handler et al. (2005) for 6* Oph (~ 8.5 Mq), where it 
seems present at a lower level. On the other hand, for V2052 
Oph we did not find such evidence, and it is hardly, if at all 
present in the MOST photometry of 7 Peg (8.5 Mq, Handler 
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et al. 2009). One might therefore speculate that the more 
massive the star, the stronger this additional variability. 

In this context it is very interesting that Blommo ct 
al. (2011) analysed CoRoT light curves of three O stars 
and also found some apparently incoherent variability, in 
all three targets. These authors suggested that it could be 
due to subsurface convection, granulation or wind variabil- 
ity. On the other hand, Balona ct al. (2011) suggested that 
the low frequencies observed in the amplitude spectra of 
Kepler B-type stars are due to many simultaneous gravity 
mode oscillations with high spherical degree. The present 
data for V986 Oph do not allow us to distinguish between 
those possibilites, and also not to argue against pulsation in 
modes of high spherical degree, because those would not gen- 
erate strong colour variability (cf. Daszyhska-Daszkiewicz 
& Pamyatnykh 2012). However, wc do point out that this 
presently unexplained variability may occur in stars with 
masses down to 9 Mq . 

Some authors have (e.g., Jerzykiewicz 1975) questioned 
the membership of V986 Oph to the class of /3 Cep stars 
due to its long variability period. Because the star rotates 
rapidly, several possibilities need to be considered. With the 
temperature and luminosity estimate for V986 Oph from the 
photometric data in Sect. 4.2 the star would have a radius of 
12 ±4Rq, which yields a rotation period of about two days 
assuming Vrot = 300 kms^^ . As the critical (break-up) ro- 
tational velocity of such massive stars is around 400kms~^ , 
the rotation frequency cannot exceed 0.7 d~^ . For the pos- 
sibility of a ~ 16 Mq, star, this upper limit increases to 
2.7 d~^ . Therefore wc rule out that the single coherent signal 
we found in the light curves of V986 Oph is due to rotation. 

Another hypothesis would be a g mode frequency ro- 
tationally split into the p/mixed mode domain. However, 
because of the large uncertainties in the stellar mass and 
effective temperature, we cannot reach a conclusion for this 
possibility and stay with the assumption that V986 Oph is 
a Cephei star. 

The frequency of the single coherent variability signal 
has changed from the first published observations, as sum- 
marized by Jerzykiewicz (1975). Up to his paper, the fre- 
quency was quoted as 3.44d~^ or somewhat higher. Later, 
PuUerton et al. (1985) gave two different periods for the dif- 
ferent seasons 1980 and 1984, the latter consistent with the 
3.29d~^ frequency dotorminod by Cuypers et al. (1989) and 
us. Therefore this frequency must have changed some time in 
the 1980's, by an amount too large to be explicable by stel- 
lar evolution. Most likely, it is just duo to a change of the 
dominant pulsation mode of the star, which has been ob- 
served in at least one other /3 Cep star before (Jerzykiewicz 
& Pigulski 1996). 

Concerning the amplitude, the published light range is 
of the order of 0.02 to 0.03 mag. This is comparable to what 
we see in our data (cf. Fig. 1). However, the amplitude of 
the main periodicity may have dropped, or the larger values 
reported in the literature, based on much smaller data sets, 
are biased by the incoherent variability. Jerzykiewicz (1975) 
already remarked on the unusually low U/B amplitude ratio 
from the viewpoint of /3 Cep pulsation, as manifested in the 
data of Hill (1967). Our observed u/v amplitude ratio is 
consistent with that, and can be best explained with an 
odd-1 mode of fairly high degree {I > 3). 

Such a mode identification is unexpected for two rea- 



sons. First, geometrical cancellation is very strong for modes 
with odd spherical degree larger than one (Daszyhska- 
Daszkiewicz et al. 2002). Second, for such a rapidly rotat- 
ing star one would naively expect a preference for Z = 2 
modes due to the distortion of the stellar shape. However, 
as demonstrated by Townsend (2003), our implicit assump- 
tion that the geometry of the pulsation mode of V986 Oph 
can be described in the form of a single spherical harmonic 
may not be correct. Also, the photometric amplitudes of 
rapidly rotating stars depend on the azimuthal order of the 
modes as well as the aspect under which they are viewed 
(Townsend 2003). 

Our most likely mode identification as Z = 3 is there- 
fore uncertain and calls for a high-resolution spectroscopic 
investigation. This would also serve to derive more reliable 
values of the stellar temperature and surface gravity than 
we have available. In particular, it would be interesting to 
confirm or reject the high stellar mass we inferred. 



6 CONCLUSIONS 

In an attempt to understand the pulsational behaviour of 

Cep stars that are more complicated than those pre- 
viously studied with asteroseismic methods (slow rotators 
with fairly large amplitudes), we have carried out an ex- 
tensive multisite campaign. Results from the photometric 
investigation were however insuSicient to perform a detailed 
asteroseismic study, as only three pulsation modes were de- 
tected for V2052 Oph, and one for V986 Oph. 

However, it is interesting that the nonradial modes 
present are of higher spherical degree than commonly found. 
We are aware of only a few cases with observationally iden- 
tified I = 3 (e.g., Briquet et al. 2009) or I = 4 modes (e.g., 
Aerts, Waelkens & De Pauw 1994), and these stars tend to 
rotate more rapidly than those dominated by modes of low 
spherical degree. 

Therefore the often-made assumption that /3 Cep pul- 
sation modes detected photometrically from the ground are 
/ ^ 2 needs to be questioned. In the context of highly sensi- 
tive space photometry this assumption is of course even more 
doubtful. Since the amplitude reduction due to geometrical 
cancellation of sufficiently high-1 pulsation modes only goes 
as ~ Z^^/^ (e.g., Daszyhska-Daszkiewicz et al. 2002), the 
gain in the number of oscillation frequencies detected is off- 
set by the larger uncertainty in mode typing. 

Photometric identifications of modes with high spher- 
ical degree become largely degenerate for even and odd 
modes with I ^ 3. This situation can be relieved by obtaining 
simultaneous spectroscopy. The amplitude ratios and phase 
shifts between the radial velocity and light curves allow some 
separation between high-1 modes (Daszyhska-Daszkiewicz & 
Pamyatnykh 2012), and of course line-profile variations of- 
fer a multitude of possibilities for identifying modes (e.g., 
Telting 2008 and references therein). 

It therefore seems that if we are to understand the in- 
terior structure of more rapidly rotating j3 Cep stars than 
those studied so far, all observationally and theoretically 
available tools need to be exploited. Photometric measure- 
ments need to be made in multiple passbands, and (simul- 
taneous) spectroscopic observations must be acquired. The 
interpretation of these data may require the inclusion of the 
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effects of rotation on a star-to-star basis because observ- 
ables can be affected to tfie extent that mode identifica- 
tions, a prerequisite for asteroseismology, may be erroneous 
(see Townsend 2003). However, to shed hght on important 
astrophysical problems, such as internal angular momentum 
transport, such concerted efforts will be worthwhile. 
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